Dilute gas of ultracold two-level atoms inside a cavity: generalized Dicke model 
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We consider a gas of ultracold two-level atoms confined in a cavity, taking into account for atomic 
center-of-mass motion and cavity mode variations. We use the generalized Dicke model, and analyze 
separately the cases of a Gaussian, and a standing wave mode shape. Owing to the interplay between 
external motional energies of the atoms and internal atomic and field energies, the phase-diagrams 
exhibit novel features not encountered in the standard Dicke model, such as the existence of first 
and second order phase transitions between normal and superradiant phases. Due to the quantum 
description of atomic motion, internal and external atomic degrees of freedom are highly correlated 
leading to modified normal and superradiant phases. 

PACS numbers: 42.50.Pq,42.50.Nn,42.50.Wk,71.10.Fd 



I. INTRODUCTION 

Progress in trapping and cooling of atomic gases [H 
made it possible to coherently couple a Bose-Einstein 
condensate to a single cavity mode These experi- 
ments pave the way to a new sub-field of AMO physics; 
many-body cavity quantum electrodynamics. In the ul- 
tracold regime, light induced mechanical effects on the 
matter waves lead to intrinsic non-linearity between the 
matter and the cavity field 0, Q- particular, the non- 
linearity renders novel quantum phase transitions (QPT) 
0, see 0. Such non-linearity, due to the quantized mo- 
tion of the atoms, is absent in the, so called, standard 
Dicke model (DM). Explicitly, the DM describes a gas 
of N non-moving two-level atoms interacting with a sin- 
gle quantized cavity mode [7]. The interplay between 
the field intensity/energy, the free atom energy, and the 
interaction energy leads to a quantum phase transition 
(DQPT) in the DM. Motivated by the novel phenomena 
arising from a quantized treatment of atomic motion, it 
is highly interesting to extend the DM to include atomic 
motion on a quantum scale, and in particular to analyze 
how it affects the nature of the DQPT. It is clear that 
such generalization of the DM results in new aspects of 
the system properties. The atomic motion is directly af- 
fected by the shape of the field induced potentials, which 
in return is determined by the system parameters and 
the field intensity. In addition to the terms contribut- 
ing to the total energy in the regular DM, in this model 
we have to take into account for the motional energy of 
the atoms. This enters in a non-trivial way since the in- 
teraction energy depends on the motional states of the 
atoms. 

The DM was first introduced in quantum optics to de- 
scribe the full collective dynamics of atoms in a high qual- 
ity cavity. The DM Hamiltonian, given in the rotating 
wave approximation (RWA), reads 



Here, the boson ladder operators a^ and a create and 
annihilate a photon of the cavity mode, the Pauli ct^- 
operators act on atom z, oj, 17 and A are the mode 
and atomic transition frequencies and effective atom-field 
coupling respectively. The DM is a well defined mathe- 
matical model for any value of its parameters. Note, 
however, that whether this model describes faithfully 
some physical situation for any choice of parameters is 
not guaranteed. For moderate A, the DM appropriately 
describes the dynamics of atoms coupled to the cavity 
field, and it has been thoroughly discussed in terms of 
collapse-revivals ^Lsqueezing [9,], entanglement ^] and 
state preparation Considerable interest, however, 

has been devoted to the normal-superradiant phase tran- 
sition [ilSlil. 

In the thermodynamic limit {N, V oo, N/V — 
const.) the system exhibits a non-zero temperature phase 
transition between the normal phase and the superradi- 
ant phase. In the normal phase, all atoms are in their 
ground state and the field in vacuum, while the super- 
radiant phase is characterized by a non-zero field and a 
macroscopic excitation of the matter. Its critical coupling 
and critical temperature are [l^ 
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As was shown in |13| . the critical coupling can be seen 
as a condition on the atomic density p — N/V . Interest- 
ingly, the DQPT is of second order nature without the 
RWA, while it is first order if the RWA has been imposed 
[m . The corrections due to the RWA to various physical 
observables have been considered [3, [l^ . 

More detailed analysis about miscellaneous aspects of 
the DQPT have been presented in numerous publica- 
tions. Especially various extensions [T7l,IT8!| as well as ap- 
proximate methods 19] concerning the DQPT has been 
outlined. Recently, C. Emary and T. Brandes applied 
the algebraic Holstein-Primakoff boson representation on 
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the DM. The method turned out to be very powerful and 
have since then been apphed frequently to the DM . 

Despite the numerous publications on the DQPT, the 
existence of this phase transition was widely discussed. 
If the two-level atoms in the DM correspond to atoms in 
a ground and excited state, and the transition is direct, 
then quantum mechanics forbids the transition. This can 
be seen either by realizing the necessity of adding the, so 
called, A^-term to the Hamiltonian, or by employing sum 
rules to bound the coefficients in the Dicke model to the 
"trivial" thermodynamical phase The argumenta- 

tion of Ref. 23] can be generalized to a quite general 
no-go theorem for the DQPT [23|, but it does not apply 
if the two level atoms in the DM correspond to atoms 
in two excited states, such as Rydberg states, or if the 
transition is not direct. 

One valuable step towards an experimental realization 
of the DQPT was taken in relation with Ref. [2l|, where 
typical experimental parameters as well as losses were in- 
cluded. These authors considered the two levels coupled 
by a non-resonant Raman transition. In such conditions 
the atom-field coupling A can be tuned more or less inde- 
pendently of the term in an effective two-level model, 
and one can reach the regime of DQPT. This paper, how- 
ever, considers a situation in which atomic motion can 
be neglected due to high temperatures, e.g. the standard 
DQPT. Alternative situation, in which the atomic mo- 
tion could be neglected, would be to consider quantum 
dots interacting with a cavity mode [l8|, (2^ . 

In this paper we extend the DM to take into account 
for atomic motion in a fully quantum mechanical descrip- 
tion. The atomic motion then introduces an additional 
degree of freedom to the problem, leading to novel ap- 
pearances of the system phase diagrams. The gas is 
assumed dilute such that atom-atom scattering can be 
neglected, and that the motion is restricted to one di- 
mension due to tight confinement in the remaining two 
directions via external trapping. Furthermore, the atoms 
are assumed trapped by the cavity field itself. Con- 
sequently, a normal-superradiant QPT is not possible, 
since for a vanishing field the atomic trapping capability 
is lost. However, we may assume a lowest bound of the 
field such that at least one bound state of the trapping 
"potential" is guaranteed. This can be achieved by an 
external pumping of the cavity, which imposes a non- 
vanishing cavity field. Our research is partly carried out 
in the adiabatic regime, motivated by the ultracold atoms 
considered and its justification is numerically verified. In 
this adiabatic regime, the problem relaxes to solving a 
1-D time-independent Schrodinger equation. In particu- 
lar we study the case of a Gaussian mode profile utiliz- 
ing this adiabatic method. The situation with a stand- 
ing wave mode profile is also considered, however using 
a full numerical rather than adiabatic approach. For a 
Gaussian profile the number of bound states is crucial 
for the thermodynamics, and we find great divergences 
from the regular DM. Among these are the existence of 
both first and second order QPT's and multiple super- 



radiant phases. The second, standing wave mode, shows 
slight similarities to the model of but the QPT is 
found to be of second order, and the PT survives for zero 
temperature and finite uj opposite to the regular DQPT. 

The paper is organized as follows. In the next sec- 
tion we present the generalized DM which includes the 
motion of the atom. The adiabatic diagonalization of 
the single particle Hamiltonians utilized for the Gaus- 
sian mode profile is introduced and the general expression 
for the partition function given. The following Sec. IIIII 
considers the situation of a Gaussian mode profile. We 
thoroughly discuss the importance of bound states. Sec- 
tion [IV] instead considers a standing wave mode profile in 
a fully numerical fashion. In the appendix, however, we 
derive analytical expressions in the regimes of tight bind- 
ing which enables us with various asymptotic properties. 
Last we conclude with a summery in Sec. |Vl 



II. GENERALIZED DICKE MODEL AND ITS 
PARTITION FUNCTION 

We consider an gas of N ultracold identical two-level 
atoms, with mass m and energy level separation hfl, in- 
teracting with a single cavity mode with frequency w. For 
a low temperature gas we include atomic center-of-mass 
motion and mode variation. In the dipole and rotating 
wave approximation, the extended DM becomes 
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Here, pi and Xi the scaled center-of-mass momentum and 
position of atom i respectively, g{x) the effective position- 
dependent atom-field coupling and V is the mode volume. 
Throughout the paper we will use scaled variables such 
that h = m = 1. The case of a single atom is given by 
the generalized Jaynes-Cummings Hamiltonian studied 
by numerous authors 25, 27]. 

In the thermodynamic limit we let y — > oo and N — > 
oo such that the atomic density is fixed; p = N/V = po- 
The partition function reads 



Tr 



(4) 



where = T and T is the scaled temperature and the 
trace is over the field and atomic degrees of freedom. It 
is convenient to perform the trace of the field in terms 
of Glauber's coherent states, a\a) = a\a). In the ther- 
modynamic limit one may replace a — > a and a* 
in the evaluation of the partition function [l^ . In other 
words; in the large atom number limit the photon lad- 
der operators, or more precisely a/^fN and o^/Vn, can 
be treated as commuting operators. Using the fact that 
atomic operators mutually commute between themselves, 
for example [ii,pj] = iSij, the partition function can be 
written 



Z = I ^ e-'^-H^ {Tr 
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where the integration is over the whole complex a-plane 
and 



h{a) 
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N 



{aa++a*a-). (6) 



In the sense of ultracold atoms as considered here, the 
kinetic energy of the atoms is assumed smaller than the 
effective potential energy. Provided that the adiabatic 
potentials do not cross, it is then legitimized to perform 
an adiabatic diagonalization of the internal states [2^ . 
In this regime, the single particle Hamiltonian relaxes to 
two decoupled adiabatic ones 
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This approximation will be imposed in the next section 
considering a Gaussian mode profile. However, in the 
proceeding section dealing with the standing wave mode 
such an approach is not advocate, since then the curve 
crossings between the adiabatic potentials break adia- 
baticity (26j . The justification of the adiabatic approxi- 
mation applied to the Gaussian mode profile will be dis- 
cussed in the end of next Section. Within this regime, the 
problem has become one of solving for the eigenvalues of 
two time-independent decoupled Schrodinger equations. 
The adiabatic Hamiltonians depend solely on the norm 
I a I and in polar coordinates the angle part can there- 
fore be integrate out. By denoting the eigenvalues E^{r) 
respectively, where r — \a\ and n is a quantum num- 
ber/numbers that can be either discrete and/or continu- 
ous, we get the adiabatic partition function 



Zad = 2 / drre 
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Without loss of generality we can choose po = 1 as it 
only scales the effective atom-field coupling. It is worth 
mentioning that the numerics deal with exponentially 
large numbers, especially for small temperatures, which 
restrict the analysis to certain ranges. 



III. TRANSVERSAL THERMODYNAMICS 



Derivation of the partition function for 
transversal motion 



A Fabry-Perot cavity has eigenmodes that are, to a 
good approximation, Gaussian in the transverse and har- 
monic in the longitudinal direction. Assuming an ex- 
ternal deep trap in the longitudinal direction and one 
transverse direction, we may consider the one dimen- 
sional problem in which the atom field coupling has a 
spatial Gaussian shape. As is well known ^25i |. and seen 



from Eq. ([7]), only atoms in the "adiabatic" internal state 
corresponding to the Hamiltonian h'^ijir) will feel an at- 
tractive potential, while the others will be scattered away 
from the cavity field. We therefore consider only a sub 
" quasi" Hilbert space containing the bound states E:^ (r) 
of 



N 



■ exp 



-2jj), (9) 



where A^, is the transverse mode width. 

In order to proceed in an analytic way, we make the 
following approximate ansatz. 



J — + — exp (-2^1 « eo + Uo^ech'iqx). (10) 

The unknown constants are determined from the con- 
ditions: (i) The two functions have the same asymptotic 
values for x — *■ ±cx3, (ii) their maximum are the same and 
(iii) they share the same FWHM. Explicitly this yields 
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The bound eigenvalues of the Hamiltonian 
are known to be 1281 



C/o(r^) sech^((7r2x) 



(12) 



-(l + 2n) + Jl + 
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Let us introduce the number of bound states, for a given 
set of parameters, as N and define the function 



N 



91 



(14) 



With this, the partition function ([5]), considering only 
bound states, becomes 



Zad — 2 



dr re 



(15) 



By the variable substitution y = r'^ /N we get 

Zad r dy e^ [-P'-v+Hai iv)]] (ig) 

Jo 
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In the thermodynamic limit, this integral is solved by the 
saddle point method [29j 

Z^,^N^ max /e^[-'5"'/+i"b^(«)l] 1 . (17) 

where Ci is a constant. Note that y has the meaning of 
scaled field intensity. 

One obstacle of the above model already mentioned 
in the introduction, is the fact that for a shallow po- 
tential well the number of bound states will vanish. In 
this limit, the cavity field can no longer serve as a trap 
for the atoms. Consequently the ground state is the one 
of zero atoms, and we cannot have a proper thermody- 
namic limit N ^ oo. We therefore add the constrain of 
a minimum of one bound state is assumed. This can be 
met experimentally by including an external driving of 
the cavity mode, so that the field is non-zero through- 
out. Thus, the "normal" phase will contain a non-zero 
cavity field which is sustained by the external pumping. 
We have numerically checked that this does not introduce 
any significant changes of the analysis. 

For the number N of bound states, we have 



(l + 2iV)<^l + ^. (18) 

Naturally, N depends on the system parameters. As N is 
an integer, a change in the system parameters may bring 
about jumps between integer numbers of N. This will 
cause discontinuities in the function gi. Letting TV = 1 
we find Uo{y) > Qy- For small fields, y — > 0, the poten- 
tial amplitude Uq vanishes and the bound states cease to 
exist. However, for small but non-zero fields, the above 
inequality may be met for large couplings A and widths 




FIG. 1: Examples of the free energy per particle fi{y) of 
Eq. (|19|l . The inserted numbers give the couplings A and the 
other dimensionless parameters are Q = lo — 1, Ax = 2 and 
T = 0.2. 



B. Numerical results 

To study Eq. (flT)) . we analyze the parameter depen- 
dence of the function 

hiy) = -PhLuy + ln[g,iy)]. (19) 

Note that fiijj) is the free energy per particle, given a 
scaled field intensity y. Let us briefly discuss character- 
istics of /i(y) before approaching it numerically. The 
first part arises from the bare field, and it is energet- 
ically favorable to have a vanishing field. The second 
part contains the atom-field interactions plus kinetic and 
potential atomic energies. The interaction energy enters 
indirectly into the atomic potential part. Increasing the 
field amplitude deepens the potential well and therefore 
lowers the energy, and it is therefore more beneficial to 
have a large field. The two terms therefore compete, and 
in particular, the location of the maximmn of fiijj) de- 
pends on the particular system parameters used. Thus, 
atomic motion, directly related to the shape and depth 
of the adiabatic potential, is a crucial ingredient for the 
QPT. If the smallest possible y maximizes the function, 
the system is said to be in a " normar phase (in quotes 
because the field is still non-zero to guarantee at least one 
bound state) , while if a non-minimal y is optimal the sys- 
tem is in a superradiant phase. In the limit of large y, 
the second term diverges as In [gi (y)] y/y while the first 
term goes as — y, and we conclude that a maximum of 
/i(y) is only obtained for a finite y. These reflections 
are numerically verified in Fig. [1] showing /i(y) for four 
different couplings A. We see that there is a critical cou- 
pling Ac for which A < Ac the system is in a "normal" 
phase and for A > Ac it is in a superradiant phase. 

In Fig. [2] we display the critical coupling Ac as func- 
tion of bj while keeping the other parameters fixed. In (a) 
we present two examples for different T and in (b) two 
examples for different A^;. For the plots, the minimum 
y is taken so that there is at least one bound state in 
the well. To the left of the curves the phase is superra- 
diant, while to the right it is "normal". The structure 
of the phase diagram is clearly different from the one of 
the regular DM in which, at zero temperature, Ac -y/w- 
For certain couplings A^ , the system is always in a super- 
radiant phase independent of ui. The location of these 
resonances are insensitive to the temperature but not to 
trap width A^,. The "sharpness" of these points makes it 
possible to have a " normal" -superadiant-" normal" QPT 
by fixing all parameters but the coupling A which is var- 
ied around As . This novel feature comes about due to the 
varying number of bound states N in the well. For say 
w w 1 and a weak coupling, the system is in the " normal" 
phase with only one bound state. As A is increased the 
system goes through a QPT into a superradiant phase. 
A closer numerical study shows that this QPT is caused 
by the sudden change of going from one to two bound 
states in the trap. As the coupling is further increased, 
the discontinuity that arose from the appearance of a 
second bound state is of less importance and the system 
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FIG. 2: The critical atom-field coupling Ac for the potential 
well (|10p as function of cj. In (a) gray curve corresponds to 
T — 0.2 and black curve to T = 0.4, while in (b) for gray 
curve Ax = 3 and for black curve Ax — 2. In both plots 
fl = 1, and in (a) Ax — 2 and in (b) T = 0.2. All parameters 
are dimensionless. 



reenters the "normal" state. Hence, the presence of a 
second bound state "forced" the system out of the "nor- 
mal" phase by inducing a sudden kink/maximum in the 
function /i (y) . When the coupling is increased even fur- 
ther, the same may happen again when a third bound 
state is introduced in the trapping potential. Eventually, 
however, the system enters a superradiant phase without 
reentering the "normal" phase. 

In order to discuss the character of the QPTs we in- 
troduce the parameter 

I^{yJi{y)= max (20) 
yo<y<oo j 

that maximizes fi{y). Thus, / is the scaled field inten- 
sity where a discontinuity in it indicates a first order QPT 
and a discontinuity in its first derivative signals a second 
order QPT. Here j/o is the lower bound which assures at 
least one bound state. The parameter / corresponding to 
the T = 0.4 curve of Fig. [2] (a) is presented in Fig. [31 We 



note that the "normal" -superradiant QPT by increasing 
A is of first order, while the superradiant-" normal" QPT 
for growing A is of second order. From the figure it is 
not clear if the second order QPT is really a QPT or a 
cross-over. Figure lU showing a slice {uj — 0.8) of Fig. [3] 
around one singularity A^ ~ 1.57, confirms that it is in- 
deed a second order QPT. Thus, contrary to the regular 
DQPT of normal-superradiance, the QPTs can be either 
of first or second order nature in this extended DM. Fig- 
ure [3] does not only reveal the type of PTs of Fig. [21 but 
also that there is a number of first order QPTs between 
various superradiant phases. 
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FIG. 3; The scaled field intensity I (|20|) as a function of A 
and ui. Here the other dimensionless parameters are f2 = 1, 
= 2 and T = 0.4. 
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FIG. 4: The scaled field intensity / of the previous Fig. [2] 
(b) as function of A when fixing the field frequency u — 0.8 
and the width Ax ~ 2. This figure clearly demonstrates the 
presence of both first and second order QPTs. 

Figures [3l and [H characterize the type of QPT in the 
X — u! phase diagram. The structure of the QPTs of the 
T ~ UJ phase diagram turns out to be equally interesting. 
One example of the T — uj phase diagram is presented in 
FigEKa). As for the X — w diagram, both first and second 
order QPT exist. FigurelHl (b), displaying the number of 
bound states, confirms that the sudden changes (first or- 
der QPTs) are due to additional bound states. The T — ui 
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phase diagrams again display several different superradi- 
ant phases. In this case, however, there are both first and 
second order PTs between the superradiant phases. In- 
terestingly, our numeric analysis indicates that the T — lo 
phase diagrams seem to be fairly independent of fi. 





FIG. 5: This figure displays in (a) the same as in Fig. |3] but 
for the T — uj plane (black curve of Fig. [Sjl , while in (b) the 
number A'^ of bound states are shown. Note in (a) the large 
number of different suerradiant phases. 



Contrary to the regular DM, in the ultracold regime 
atomic motion plays an important role for the system 
characteristics. The confining potential is determined for 
a given field intensity y, and consequently, the atomic 
density Pat{x) per particle depends as well on y. In other 
words, apart from changes in the field intensity and in 
the atomic inversion, the phase transition will as wU be 
manifested in the atomic density. A natural consequence 
of the coupled system considered here is that the mo- 
tional state of the atoms are entangled with the internal 
state. Thus, the phases (normal and superradiant) arc 
intrinsically different from the ones of the regular DM. 
The regular DQPT derives from a competition between 
the free field energy and the interaction atom-field en- 
ergy. While the free field is minimized by vacuum, the 
interaction energy decreases with an increasing field in- 



tensity. In the present model does the kinetic energy 
contribute to the total energy. It is thus an interplay be- 
tween three terms; free field, atom-field interaction (con- 
taining the bare internal atomic energies), and motional 
energies. This is not so evident from the free energy per 
particle where the kinetic energy is hidden in the 

second term. To illuminate the importance of the atomic 
motion we study the atomic inversion W defined as the 
probability for a single atom to be in its excited state 
|e) minus the probability for it to be in the ground state 
\g). Once the adiabatic approximation has been imposed 
we are left with a single internal state; the lower adia- 
batic state \—)ad = single) -l-cos^jg). In particular, the 
angle tan 20 = 2g{x)yfpQa/Q.y/N depends on the spatial 
coordinate and the inversion becomes 

ly = Tr [a,pat{x)] ^ J dx (sin" 9 - cos" 9) pat{x). (21) 

The above equation clarifies that the atomic motion en- 
ters the problem in a non-trivial way. It also shows 
how the internal atomic properties are taken care of even 
though in the adiabatic approximation the system prop- 
erties derives from a single internal state \—)ad- We have 
numerically verified the appearance of the PT in terms 
sudden changes in Pat(x) (first order PT) or dpat{x) / dx 
(second order PT). 

In deriving the phase diagrams, tight confinement of 
the atoms in two directions has been assumed. If only the 
longitudinal motion is frozen out, one regains an effective 
two dimensional problem whose eigenvalues are obtained 
from the Schrodinger equation with potential V{x, y) — 

— Aexp {^— ^^i y As in the one dimensional situation 

studied in this section, V{x,y) possesses a finite number 
of bound states and one would expect very similar phase 
diagrams for this two dimensional case as for the one 
dimensional model. 



C. Validity of the adiabatic approximation 

We conclude this section by analyzing the adiabatic ap- 
proximation. By a simple rotation, the amplitudes a ap- 
pearing in the single atom Hamiltonian ([6]) can be taken 
real. For real a, the two last terms of h{a) is readily 
diagonalized by the unitary transformation [26| 



U 



cos 9 sin 9 
— sin 9 cos 9 



(22) 



where the angle 9 was given right above Eq. (pij) . Mo- 
mentum transforms as UpW = p — {a^ + '^^) dS, d9 = 
d9/dx. Due to the spatial dependence of 9 = 9{x), 
the transformed Hamiltonian is non-diagonal; h{a) = 
Uh{a)W = had{\a\) + hcor{a). Here, hadi\a\) is the 
adiabatic Hamiltonian ([7]) and hcorict) contains the non- 
adiabatic corrections. Explicitly one finds [1^ 



id9y 

-2i {89) p 



2i {89) p + 8^9 
{89)^ 



(23) 
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The kinetic energy is smaller or of the same order as 
maxx V^{x, jap), which provides a measure of /icor (a) in 
comparison to the adiabatic Hamiltonian /iad( |a|)- For 
typical parameters, a; = 17 = 1, A^; = 3, and y — 2, 
the terms of hcorict) are at least one order of magnitude 
smaller than the terms oihadi\c(\), which justifies the use 
of the adiabatic approximation. 



IV. LONGITUDINAL THERMODYNAMICS 

In the previous section we studied how the DQPT was 
modified due to motion of the atoms in a finite potential 
well, assuming the atomic motion to be frozen out in 
the longitudinal and one transverse direction. Here we 
instead assume the atoms to move freely along the center 
axis of the Fabry-Perot cavity while tightly bound in the 
transverse directions. 

The corresponding single atom Hamiltonian ([6]) reads 



f nn^^ Acos(mx)^ ^ 

n{a) = 1 a + n ^ [acr^ + a a 

2m 2 ^/n ^ 

(24) 

where ^ is the scaled photon wave number which will be 
set to unity hereafter, /i = 1. This Hamiltonian, with 
the field still quantized, has been considered in several 
papers, see for example [131 . Normally L ^ 2tt, where 
L is the cavity length, so neglecting boundary effects is 
not a crude approximation [27| . The Hamiltonian is of 
the form of a generalized Mathiue equation, and hence 
the spectrum E^{k) is described by a band index v and 
a quasi momentum k extending over the first Brillouin 
zone. Due to the internal two-level structure of the atom, 
the Brillouin zone is twice the size of what is imposed by 
the periodicity of the mode Clearly, E,^{k) depends 
on the field amplitude a. The corresponding eigenfunc- 
tions are written ^k,uix) = 4''ktix)\e) + (j^'kllg)- For a 
constant coupling g{x) = Aq, these Block functions are 
simple plane waves giving a constant energy shift in- 
dependent of system parameters such as A, w, and Q. 
For a standing wave mode coupling on the other hand, 
the Bloch functions cannot be decoupled from the inter- 
nal atomic states, and consequently the atomic motion 
will affect the structure of the phase diagrams as will be 
demonstrated below. 

In the previous section we assumed the adiabatic 
regime and diagonalized the Hamiltonian in its internal 
degrees of freedom. In this case, the adiabatic potentials 
V^{x) cross and adiabaticity breaks down in the range 
where the QPTs occur. Fortunately, the Hamiltonian is 
easily diagonalized numerically by truncation the dimen- 
sion of the Hamiltonian matrix. We present, however, 
asymptotic analytical results in the Appendix which re- 
lies on the adiabatic approximation. These analytical 
results enable us to extract the limiting situation of large 
field amplitudes. Furthermore, as a numerical diagonal- 
ization directly renders several of the Bloch bands we do 
not restrict the analysis to just the lowest one. However, 



it turns out that for most of the presented examples only 
the lowest band contribute to the dynamics due to the 
low temperatures considered. Exceptions are in the plots 
of the critical temperature where we indeed go to rather 
high temperatures and the excited bands become impor- 
tant. 

The partition function is written like in the previous 
section as 



Z : 



7V-S= max <!e^^^(^) 



_/V 0<j/<oo 



(25) 



where C2 is a constant and the free energy per particle 
f2iy)^-phLoy + ln[g2iy)] (26) 

with 



92 



00 „ + ! 



(27) 



Here, as above, y = \a\'^ /N represent the scaled field 
intensity. Shown in the Appendix, the second part of 
72(2/) scales asymptotically as ~ ^ for large y. Thus, 
a maximum of the free energy can only be obtained for 
finite or zero field intensities y. 

As in Sec. lIIIi we derive the critical atom- field coupling 
Ac and temperature T^. In Fig. [5] we show the results of 
how the critical coupling depends on uj (a) and (b) 
for different temperatures. The critical temperature as 
function of u and 17 is displayed in Fig. [71 where the in- 
serted numbers indicate the values of the coupling A. In 
both cases, the critical quantities show clear differences 
compared to the ones of the regular DM, ([2]). The crit- 
ical coupling scales as Ac ~ for fixed O just like in 
the regular DM. However, the H-dependence is not pos- 
sessing the same structure as dH), and especially for small 
values on 17 the critical coupling Ac is non-zero. This was 
also found in a nearest-neighbor coupling model studied 
in [3. 

In the regular DM, for a fixed A and 17 the critical 
temperature diverges for small w and goes to zero for 
large u. Here we note that for high temperatures the 
regular behavior is regained, while for low temperatures, 
and in particular zero temperature, a QPT takes place 
for finite uj. Fixing oj instead and vary 17 we get even 
more surprising results. In the regular DM, there is an 
upper temperature for which the QPT is lost and at zero 
temperature a QPT occurs for finite 17. In our model, 
a similar phase-diagram is obtained for a range of pa- 
rameters OJ and A, but there also exist parameter regimes 
where no QPT occurs for zero temperature. 

Like in the previous section the nature of the PTs is 
studied by introducing the scaled field intensity / max- 
imizing /2(y)- It is found that the QPT's are of second 
order character in all cases. 



8 



T=3 / 












/ 






■ / / 






1 / 
1 / / 







0.2 0.4 0.6 0.8 

CO 




FIG. 6: The dimensionless critical atom-field coupling Ac as 
function of uj (a) and Q (b). In (a) = I and in (b) cu = 1. 
The numbers to each curve display the dimensionless tem- 
perature T. Note in particular that for $7^0, the critical 
coupling Ac 7^ 0. 
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FIG. 7: The dimensionless critical temperature Tc as func- 
tion of the system parameters. The inserted numbers give 
the value of A. We note that a QPT is possible at zero tem- 
perature and finite uj, and moreover, that for zero Q. the QPT 
may vanish for small temperatures. In (&) Q, = 1 while in (b) 



V. CONCLUSIONS 

In this work we have studied a new regime in the DM. 
The atoms are assumed trapped by the cavity field itself 
and ultracold such that their center-of-mass kinetic en- 
ergy is of the order of the atom-field interaction. This 
calls for a full quantum mechanical treatment of the 
atomic motion and at the same time take into account 
for spatial mode variations. The analysis is motivated 
from our earlier findings, where we demonstrated that 
atomic motion greatly affects the system dynamics in 
many-body cavity QED systems [3| ■ Expectedly, we have 
shown that this is also true for the DM. The analysis is 
restricted to considering one dimensional problems, and 
both the case of a Gaussian and a standing wave mode 
profile were treated. However, we motivated that simi- 
lar phenomena are expected also for higher dimensional 
situations. In particular, we made evident that the vary- 
ing number of bound states in the potential formed by 



a Gaussian mode profile induces novel first order QPTs. 
Additionally, we found that for certain couplings As, the 
system is superradiant for any field frequency to. More- 
over, great differences with the regular DM was also en- 
countered for a standing wave mode profile. 

APPENDIX A: TIGHT BINDING APPROACH 

Here, we analytically consider the large field asymp- 
totic expressions for the energy per particle f2{y) in the 
case of a standing wave mode profile. In this regime one 
may utilize the tight binding approximation [soj to derive 
the spectrum. We further assume the adiabatic approxi- 
mation to be valid and that we can restrict the analysis 
to the lowest energy band. The adiabatic potentials are 
given by 



V^^ix) =±h^— + A2cos2(x)2/, (Al) 
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where y = \a\'^ /N is as before the scaled field inten- 
sity. A convenient base for writing down the periodic 
Hamiltonian in matrix form is to use the Wannier states 

'^fip^) = {x\j)w,± [13] ■ The function wf{x) is localized 

' -I- 
in the jth "well" of the potentials V^d(^)- tight 



binding approximation we assume w,±{'i'\h^ii\j) 







unless i = joTi = j±l. Within the validity regime of 
this approximation we may as well replace the Wannier 
functions by Gaussian functions [1]. The widths of the 
Gaussians are given by approximating the potential wells 
by harmonic oscillators, giving 



wf{x) W Wq{x — Xj) 



1 



(A2) 



The potential integrals of (jA4p are not analytically solv- 
able for the given potentials (|Aip . Instead we make the 
same kind of approximation as in section [TTTl 



and identify 



V^^{x) « ±A±B cos'^ix), 



2' 



(A7) 



(A8) 



where 



Within this regime we find 



(x) 



dx^ 



(A3) 



and Xj is the position of the jth potential well. To avoid 
un-physical contributions from the non-orthogonality of 
the Gaussians we impose / dxwG(x—Xj)wG{x—Xi) — Sij. 
We further introduce the matrix elements 



My)-' 



4(x- 



Xj+i) 



dxW(j {x-Xj)V^^{x, \a\^)W(j{x~Xj+i), 



(A4) 

where we only consider i = 0, 1. We note that the Wan- 
nier functions are directly related to the depth of the 
corresponding potential and therefore their width a will 
also depend on y. This explains the field intensity depen- 
dence of Ei{y). Another important observation is that 
WQ{x — Xj) are localized where |V^'rf(a;)| are close either to 
its maximum or its minimum, resulting in different cou- 
pling elements J^{y), indicated by the zb-superscript. In 
this notation we get the lowest band tight binding energy 



Et{k) = Eo{y) + 4{y) 



[Ei{y) + jt{y)]2cos{k). 



(A5) 



The part in front of the cosine function is strictly nega- 
tive resulting in that the ground state energy is given by 
fc = 0. The kinetic energy integrals of (|A4p are readily 
solvable, and one finds 



i^i(j/)--^exp(--^)(2a^+.^). 



(A6) 



2 Ay/Wia^ 



Ji (2/) 3^e 



(A9) 



We emphasize that the width depends on the field 
intensity y; 



1 

2B' 



(AlO) 



The applied approximations are only reliable for z < 1 
^4], and it turns out that the QPTs occur beyond these 
approximations. Nonetheless, we may find the asymp- 
totics for the free energy /2(y). In the large y limit we 
find that In [32(2/)] ~ \/y- Consequently, the field in- 
tensity / will always be finite, regardless of parameter 
choices. We have verified numerically the y square-root 
dependence of In [52(2/)] for large intensities. 
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